Implementation of a Blowing Boundary 
Condition in the LAURA Code 
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Preliminary steps toward modeling a coupled ablation problem using a finite-volume 
Navier-Stokes code (LAURA) are presented in this paper. Implementation of a surface 
boundary condition with mass transfer (blowing) is described followed by verification and 
validation through comparisons with analytic results and experimental data. Application of 
the code to a carbon-nosetip ablation problem is demonstrated and the results are compared 
with previously published data. It is concluded that the code and coupled procedure are 
suitable to support further ablation analyses and studies. 
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Nomenclature 

= blowing parameter 
= species mass fraction 
= Blasius function 
= static enthalpy 
= heat transfer coefficient 
= total enthalpy 
= recovery enthalpy 
= reference length 
= mass flow rate 
= pressure 
= heat transfer rate 
= Reynolds number 
= sphere nose radius 
= surface distance 
= Stanton number 
= time 
= temperature 

= streamwise and normal velocity 
= coordinate directions 
= correction factor 
= Blasius coordinate 
= density 


Subscripts 
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species number 
boundary layer edge condition 
stagnation condition 
reference condition 
streamwise direction 
wall condition 

condition with zero mass transfer 

condition at volume face one-half cell above wall 

condition at volume face one-half cell below wall 
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I. Introduction 

T he Crew Exploration Vehicle (CEV) is being currently designed to support manned missions returning from 
low Earth orbit and the moon. At lunar return velocities, all currently available thermal protection systems 
(TPS) will rely on an ablative heatshield design to absorb the high energies associated with Earth reentry. 
The modeling and design of these heatshields is traditionally done using an uncoupled approach where the heating 
environment is computed in the absence of ablation and is later modified by a blowing correction to determine the 
TPS material response. However, some efforts at performing a coupled analysis have been reported (Refs. 1-3) 
where computational solutions of the flowfield and the ablating material were performed in an iterative fashion. 
While much success has been achieved using an uncoupled analysis for design, it is of interest to obtain higher 
fidelity (coupled) analyses in order to better assess the system and quantify design margins. This paper will present 
steps toward those goals by documenting the implementation of a mass-transfer surface boundary condition 
(blowing) into a state-of-the-art code (LAURA, Ref. 4) that can be used in future coupled analyses. Similar progress 
has been reported by Martinelli, et. al. (Ref. 5) for another computational fluid dynamics (CFD) code. 
Modifications to the LAURA code needed to model an ablating surface are described in this paper followed by 
verification and validation of the results. While mass injection and ablation processes are often turbulent, only 
laminar flow conditions are considered in this paper. A process for coupling the LAURA code in a surface ablation 
and recession analysis will be described and results for a carbon nosetip material will be demonstrated. 


II. Implementation 

The LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) code (Ref. 4) uses a finite-volume 
shock-capturing approach to solve steady viscous and inviscid flow problems. The algorithm incorporates point- 
implicit or line implicit relaxation schemes to obtain solutions efficiently on multi-processor and massively parallel 
computers. The code has been successfully applied to a wide range of hypersonic vehicles and flight conditions 
during its evolution over the past 15 years. LAURA includes models for perfect gas, equilibrium air, and thermal 
and chemical nonequilibrium air in addition to models for CF 4 and a Mars atmosphere. An important feature of the 
code is the ability to perform one-dimensional grid adaption in parallel with the solution to resolve high gradients in 
the boundary layer and across a bow shock. For the applications in this work, the thin-layer form of the Navier- 
Stokes equations were solved and an eigenvalue limiter of 0.3 was used in the code. The computational grids in 
every case were adapted such that the cell Reynolds number at the wall was 0(1) and solution convergence was 
obtained when changes in surface heating were less than 0.1% after 5000 iterations. 


The boundary condition for modeling gas injection at a solid surface is fairly straightforward and was treated in 
an explicit fashion in the LAURA code. Figure 1 illustrates a computational cell centered on a wall boundary with 
flow conditions one-half a cell above the surface (point 2) and one-half cell below (point 1). This lower cell is 
treated as a a pseudo-cell for the purposes here. At the wall, the mass flow rate (p K v u ,). temperature (T w ), and species 
mass fractions (c,j are assumed to be known from a material response code and by assuming chemical equilibrium at 
the surface. 
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Figure 1. Cell volume for treatment of wall blowing boundary condition 
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Assuming no losses through the vertical (side) faces, the conservation of momentum in this cell can be written 


Pi+pi^i = Pw +py w 

(1) 

with the massflow at the wall 


m = p v 

rw w 

(2) 

and the pressure and density related by thermodynamics 


P W = Pn(PnJn) 

(3) 


This system of three equations in three unknowns can be solved to yield the density, pressure, and velocity at the 
surface. In this explicit boundary condition treatment, a solution for the surface quantities is performed at each 
iteration step in LAURA. During development of the boundary condition a simple test was devised using a sphere 
geometry in Mach 6 air with helium injection at the surface as shown in Figure 2. 


Contours of He mass fraction 



Figure 2. Test case for development of blowing boundary condition 

In the test case, a region of mass injection was defined beginning at the stagnation point and ending a short distance 
downstream. A range of mass flow rates were simulated to test the implementation of the blowing boundary 
condition and study solution convergence. The figure illustrates the flowfield computed for a case with a moderate 
blowing rate where the color contours represent the different gases — helium (red) and air (blue). Streamlines from 
the surface injection and oncoming flow illustrate the computed velocity field. These qualitative results appear as 
expected and variation of mass flow rates proved the boundary condition scheme to be robust within LAURA for the 
cases considered. 


III. Results 

The following sections present results that were achieved with the code modifications in order to verify and 
validate (V&V) the blowing boundary condition implementation. This V&V process is expected to be a continuous 
task as new conditions are analyzed and new results are identified that can be used for comparison. 

A. Verification 

Initial verification of the blowing boundary condition in LAURA was performed by using the code to compute 
the theoretical Blasius velocity profile for flow over a flate plate with injection. In the Blasius approximation, the 
equation for flow over a flate plate is given by 

/"+#' = o < 3 4 > 
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where the boundary condition for mass flow through the surface, in terms of the Blasius function, is 

P W V W _ /( 0 ) 

Poo v oo 27Re7 (5) 

A rectangular computational domain was constructed to compute the Blasius flat plate problem with LAURA. The 
domain extended 10 m in the vertical direction and 10 m from either end of a 1 m plate in the upstream and 
downstream directions. As constructed, the CFD problem included the effects of flow moving onto and off the flat 
plate at the leading and trailing edges, respectively. Figure 3(a) illustrates those effects on the velocity profiles near 
the plate edges (x=0 m and x=1.0 m) for the case with no blowing through the surface. As shown in the figure, the 
velocity profiles near the center of the plate (0.5 < x < 0.75 m) collapse to the Blasius profile when plotted in the 
similarity variables. A similar plot is shown in Figure 3(b) for additional cases where blowing into the boundary 
layer is introduced at the surface. The velocity profiles in Figure 3(b) correspond to cases where the Blasius 
function at the surface, /(0), equals 0.0, -0.4, and -0.8. Flere, the negative values for/(0) denotes flow into the 
boundary layer (blowing); no solutions for positive values of f(0) (i.e. for suction) were computed. The data used 
for comparison in Figure 3 were obtained from the tabulation of Blasius profiles reported by Emmons and Leigh 
(Ref. 6). Good agreement between the LAURA computed profiles and the Blasius data was found. 



Figure 3. Comparison of predicted velocities with Blasius profiles 


An early approach to the analysis of laminar boundary layers with mass transfer was the use of film theory such as 
reported by Mickley, et.al (Ref. 7). Following the thin-film approximation, it was shown that variation of the 
transfer coefficients (friction, heating, and diffusion) can be represented by a simple exponential function 


6 = 



(6) 


where 6 is a ratio of the coefficient with and without mass transfer and (f) is a function of the mass transfer rate. 
Derivations of the same function have been presented by several authors (Refs. 8-10) and the resulting function for 
the reduction of heat transfer in the presence of blowing has been widely used. For heat transfer, the function in 
Eq. 6 is usually written in the form 

St B 0 

= (7) 

where the Stanton number is defined by e 1 


St = 


q 

Pe U e( H r~K) 


( 8 ) 
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and the blowing parameter by 


( 9 ) 


B 


0 


Pe U e St 0 


The subscript 0 in the preceding equations denotes a quantity taken at a zero mass transfer condition. The 
assumptions associated with Eq. 4 are that the flow has Couette-like properties so that gradients in the streamwise 
direction are much smaller than across the layer. In addition, the derivation assumes that the mass flux introduced at 
the surface is constant across the layer and that the layer thickness is unaffected. There is no explicit requirement 
(or modeling) for the flow to be laminar or turbulent in the above equations. A slightly modified version of Eq. 4 
has the form „ „ „ „ 

St _ 2AB 0 

~St~ 0 ~ e 2AB ° - 1 ( 10 ) 


or, equivalently 


St ln(l + 2 AB) 
~St^~ 2A B 


( 11 ) 


where 



Pe U e St 


( 12 ) 


Introduction of the A term in Eq. 10 was reported (Ref. 11) to improve the agreement with turbulent data when 
A= 0.4 while reducing to the original formula for laminar flow when T=0.5. Equation 11 is the basis for correcting 
the surface heat transfer in the CMA (Ref. 11) and FIAT (Ref. 12) material response codes when an unblown value 
of heating ( St 0 ) is supplied. 


To provide further verification of the boundary condition implementation in the LAURA code, a model problem 
was solved for a Couette flow using the assumptions described above. It was expected that the reduction in surface 
heat transfer computed with LAURA should match the values predicted by Eqs. 7-9. For the CFD solution, a 
rectangular domain was constructed with a planar surface representing the wall and freestream conditions prescribed 
at the outer (edge) location. The upstream and downstream faces of the domain were set to match each other since 
the flow was assumed constant in the streamwise direction. Perfect gas flow solutions were obtained for a range of 
edge conditions that included Mach numbers from 0.15 to 3.0, unit Reynolds numbers from lx 10 4 to lx 10 5 m , and 
wall temperature ratios ( TJT„) between 0.2 to 0.65. Results for the computed variation in heat transfer for blowing 
parameters ( B ) between 0 and 3 are shown in Figure 4. These computed values from LAURA are compared with 
the analytic result from Eq. 7 and, while displaying some disparity, clearly show that the level and trend are very 
closely predicted. 



Figure 4. Stanton number reduction as a function of blowing parameter 
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B. Validation 


While a substantial amount of theoretical work exists in the literature concerning the effects of blowing on heating 
rate, there are a limited number of experimental tests reported that can be used for code validation. In the present 
work, comparisons have been made with the published results of Feldhuhn (Ref. 13), Kaattari (Ref. 14), and Marvin 
(Ref. 15). Additional sources of data from Libby (Ref. 16) and Laganelli (Ref. 17) were identified during the course 
of this work but those data were dominated by turbulent flow and comparisons are not presented in this paper. 

Feldhuhn measured heat transfer distributions on a porous, spherically blunted cone in Mach 5 air over a 
substantial range of Reynolds numbers. The 5-deg sphere-cone model was constructed with three interior manifolds 
such that the mass transfer rates of the injected nitrogen could be varied independently in the stagnation region, the 
sphere shoulder, and the conical flank. A matrix of tests were performed where the blowing varied between zero 
and seven times the nominal to yield injection rates up to 2 percent of the freestream massflow. At the higher mass 
flow rates and Reynolds numbers, the test data clearly showed transition to turbulence at the location of mass 
injection or downstream on the conical flank. A small number of wind tunnel conditions were tested where the flow 
was believed to remain laminar over the sphere cone nose and the data from those tests were used for validation in 
the present work. 

Four test conditions were considered for comparison (Cases 15, 22, 23, and 141 in Ref. 13) with each case being 
at Mach 5 and 50 atm total pressure which yielded the lowest freestream Reynolds number (3.8 x 10 6 ft" 1 ). Mass 
flow in the stagnation region and sphere shoulder was set at zero, once, or twice the nominal value (0.6 percent of 
freestream). These varying injection rates were defined in LAURA for the CFD computation and a comparison of 
the measured and predicted heating for these cases is shown in Figure 5(a) and 5(b). The results computed in this 
work have been cast in terms of heat transfer coefficient and plotted on the same scale and axes as the data was 
presented in Ref. 13. 




b) Blowing over the entire surface 


Figure 5. Comparison of heat transfer coefficients on a sphere for different levels of mass transfer 


The surface blowing in Fig. 5(a) corresponds to cases where there was zero mass transfer in the stagnation region 
(Region 1) prior to a surface distance of s/R n =0.52. Downstream of that location, the mass transfer was set at 1- 
times and 2-times the nominal rate which resulted in reduced surface heating over the sphere shoulder (Region 2) as 
shown in the figure. In Fig 5(b), a nominal blowing rate was introduced in the stagnation region while the 
downstream mass transfer rates were approximately the same as in Fig. 5(a). With blowing over the whole surface 
(Fig. 5b), the heating is lowered by ~30% at the stagnation point with additional reduction over the shoulder 
depending on the blowing rate. In both figures, there is good agreement between the measured heating and the 
values computed with LAURA. Some discrepancy with the data is noted around s//R n =1.2 but the measurement 
appears spurious at that station. 


The tests by Kaattari (Ref. 14) consisted of a hemisphere geometry and a blunted spherical segment 
representative of a capsule heatshield and so are of considerable interest in relation to CEV. Additional tests of a 


6 

American Institute of Aeronautics and Astronautics 





blunt conical body were reported in Ref. (14) but those data were not studied in this work. Figure 6 presents the 
experimental setup for the two models considered herein. Each model consisted of a permeable header in the shape 
of the hemisphere or spherical segment which was attached to a pressurized chamber in the body of the model. 
Mass flow rates through the permeable surface were calibrated during pretest based on the tunnel total pressure, 
surface pressure (predicted or measured), and flow conductance of the porous surface. A result of the calibration 
was that the mass flow varied along the surface with a lowering of injection rates in the stagnation region of the 
hemisphere model. The tabulated mass flow rates from the calibration and reported in Ref. 14 were matched exactly 
in the CFD calculations. 


CALORIMETER 

DETAIL 



Figure 6. Blowing test model for hemisphere and spherical segment from Ref. 14 

Predictions of the surface heat transfer using LAURA for the two shapes studied are compared to the measured 
heating rates at different levels of blowing in Figs. 7(a) and 7(b). Freestream conditions for both cases were at Mach 
7.32 and a total pressure of 552 N/cm 2 . At the lowest level of blowing, the prediction and measurement are in good 
agreement on the hemisphere but have increasing disparity in the stagnation region as the mass flow is increased. In 
the case of the spherical segment (Fig. 7b), the agreement between measurement and prediction is not good in the 
stagnation region but improves substantially in areas away from stagnation region for both levels of mass transfer. 
In these comparisons, the data measurements indicate that the reduction in heating at the stagnation point is 
substantially less sensitive to the mass transfer rate than was predicted by LAURA. The reason for disagreement 
between some of the data and the predictions is not completely understood, however, it was noted in Ref. 14, that 
the data from Feldhuhn (Ref. 13) suggests a similar level of insensitivity when the flow in the injection region was 
turbulent. Any such turbulent flow effects not modeled in the computations coidd explain the discrepancies. 
Overall, the comparisons between prediction and measurement were inconclusive for these data. 
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Figure 7. Comparison of heat transfer ratio with mass transfer for data of Ref. 14 
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The data published by Marvin and Akin (Ref. 15) offers different flow physics for validation of the blowing 
model in LAURA. In that experimental work, a 5-deg cone was fitted with a impermeable nosetip for the first 3.75 
in of vehicle length followed by a porous conical section. Heat transfer measurements were made using air, argon, 
and helium injection along the cone surface although only air injection was considered in this work. The model was 
fitted with both sharp and blunt noses while freestream Reynolds number and injection rate were varied to 
investigate the effects of those parameters on boundary layer transition. Since the experiment was designed to 
promote transition on the conical flank, an adequate extent of laminar flow existed for most cases to make 
comparisons between measurement and LAURA predictions. The metallic experimental model had a uniformly 
porous surface that was calibrated so that mass flow rates along the vehicle length were accurately known. The 
calibration yielded a mass transfer distribution with ~20% variation about the mean level. The reported distribution 
was input to LAURA as the surface boundary condition and calculations were done with four mean levels: 0%, 
0.04%, 0.08%, and 0.12% of the freestream mass flow rate. The freestream Mach number was 7.4, total temperature 
was 600 psi, and Reynolds number based on length was 4.7 x 10 6 for the cases computed. A comparison of the 
predicted heating rates with the experimental data is presented in Fig. 8 for the different blowing levels. In this 
figure, the heating rate is normalized by the experimental value of heating at the beginning of the porous cone 
section and plotted in the same coordinates and that were used in Ref. 15, however, some of the turbulent data that 
was measured is not included here for clarity. The experimental data in this figure clearly exhibits a transition to 
turbulent levels along the cone flank with the exact location dependent on the level of blowing. The transition is 
most forward at the highest levels of injection as would be expected. Excellent agreement between the measurement 
and computation was found for these data and served to further validate the boundary condition implementation in 
LAURA. 



Axial Length, s/L 

Figure 8. Heating rate ratio for sharp cone with air injection 
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IV. Application 

Since the aim of implementing the blowing boundary condition into 
LAURA is to enable coupled ablation analysis it was of interest to 
demonstrate this capability for a test problem. Residts have been 
published (Ref. 2) using a suite of codes to analyze the time- 
dependent ablation of a 7-deg sphere -cone reentry 

vehicle (RV) with a 1.5 inch radius carbon nosetip. The simplified 
geometry, single carbon ablation element, and extensive published 
results offer a good benchmark for a demonstration problem. Kuntz, 
et.al. (Ref. 2) used a 2-dimensional in-depth material response code 
(COYOTE) in their work and allowed for shape change of the nosetip 
during the analysis. In the present work, the nose shape was allowed 
to vary in response to the ablation, however, a one-dimensional in- 
depth response code (FIAT, Ref. 12) was employed and so the 
material response was expected to vary somewhat from previous 
results. Another difference in the present work in comparison to 
Ref. 2 was that a simple marching of the coupled problem was 
employed herein, whereas a predictor-corrector scheme was employed 
by the authors in Ref. 2. In their approach, the aerothermodynamic 
environment for each time step was computed by a CFD code and 
then corrected by iteration at each time step. It was reported (Ref. 2) 
that this approach was necessary for stability of the coupled analysis 
due to the high rate of ablation (and shape change) at later times along 
the RV ballistic trajectory. In the present work, the 
aerothermodynamic environment was computed at ~1 s time steps and 
held fixed while the material response was analyzed over that period. 
The whole analyses was ended before the period of massive ablation 
was reached to avoid stability problems. The overall procedure 
employed for the coupled analysis is diagrammed in Fig. 9 where the 
various codes employed are named in the boxes and the arrows 
between boxes show the dataflow from code to code. This procedure 
was essentially the same as reported in other coupled analyses (Refs. 
1,2). The ACE code (Ref. 18) in Fig. 9 was used to compute the 
equilibrium species at the surface once the pressure and temperature 
were known from the tlowfield (LAURA) and material response 
(FIAT) codes. Shape change of the vehicle geometry was determined 
from the recession computed in FIAT and regridding for the next 
flowfield solution was done using a code (MVGRD) written for this 
test problem. Figure 10 presents the original surface and final ablated 
shape and associated flowfield grids after a total of ten time steps 
were analyzed following the methodology described. 



Figure 9. Coupling procedure for 
ablation anlaysis 
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1.5 i— 


t = 0 sec 


Initial body shape and flowfield grid 



Figure 10. Initial and final (ablated) surface and flowfield grids 

Figures 11(a) and (b) compares the predicted stagnation point recession and wall temperature in comparison to the 
results from Ref. 2. Good agreement with those published data was found thereby demonstrating that a coupled 
analysis capability using the LAURA code to provide the aerothermodynamic environment with ablation mass 
transfer at the surface. 



Figure 11. Comparison of coupled ablation predictions with results of Ref. 2 
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V. Conclusion 

The implementation of a boundary condition with mass transfer (blowing) in the LAURA code proved to be a 
straightforward task and a range of mass transfer rates have been modeled in test cases with air being the primary 
injectant. Comparisons of the predicted velocity profiles with tabulated Blasius data showed good agreement as 
did a comparison of the reduction in Stanton number due to mass transfer from an analytic Couette flow solution. 
The capability to predict heat transfer in the presence of blowing was validated by comparisons with blunt nose 
and sharp cone data. Solutions with finite-rate chemistry in the flowfield and injection of carbon at the surface 
were also obtained in a coupled analysis of ablation and surface recession. Comparison of those results with 
previously published data for wall temperature and nosetip recession showed good agreement and demonstrated 
new capability. Computer calculations for the range of cases considered have shown that the CFD solutions 
behaved robustly and convergence was easily achieved. Computational performance with more complex surface 
chemistry (i.e. more and varied injectant species) remains to be measured but it was determined that the code and 
coupled ablation procedure is suitable for further applications. 
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